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ABSTRACT 

We examine the biases induced on cosmological parameters when the presence of secondary 
anisotropies is not taken into account in Cosmic Microwave Background analyses. We first 
develop an exact analytical expression for computing the biases on parameters when any 
additive signal is neglected in the analysis. We then apply it in the context of the forthcom- 
ing Planck experiment. For illustration, we investigate the effect of the sole residual thermal 
Sunyaev-ZeFdovich signal that remains after cluster extraction. We find in particular that 
analyses neglecting the presence of this contribution introduce on the cosmological param- 
eters n s and t biases, at least ~ 6.5 and 2.9 times their one cr confidence intervals. The fib 
parameter is also biased to a lesser extent. 
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1 INTRODUCTION 

Future Cosmic Microwave Background (CMB) experiments, which 
are designed to be cosmic variance limited, will allow us to deter- 
mine the cosmological parameters with a relative precision of the 
order of, or better than, one percent. It will be made possible in 
particular through the measurement of CMB temperature and po- 
larisation anisotropies with unequalled sensitivities, exquisite an- 
gular resolution and optimal frequency coverage. In this context, 
additional contributions to the signal (galaxies, point sources, sec- 
ondaries arising from the interaction of CMB photons with matter 
after decoupling, etc.) can no more be neglected. More specifically, 
a precise quantification of the biases on the parameters and of their 
sources is now needed. 

The study of biases in cosmology is receiving growing 
attenti on. In the context of weak lensing, lAmara & Refregien 
d2008h have derived a method based on a Fisher matrix type 
analysis for quantifying system atic biases. In CMB analyses, 
Miller, Shimon & Keating (2008) examined the biases introduced 
by beam systematics for five upcoming experiments that will 
measure the B-mode polarisation (Planck, PolarBeaR, Spider, 
Q/U Imaging Experiment (QUIET)+Clover and CMBPol). Simi- 
larly, using Planck characteristics, previous studies have estimated 
the biases induce d by the contribution from patchy reionisation 
dZahn et alj|2005h . They concluded that the biases, depending on 
the model of reionisation, can be as high as a fe w in units of the 
one sigma error. More recentlv. lSerra e t al. (2008) have focused on 
the contribution from clustered IR point sources and its effect on 
the cosmological parameters. Those two studies used quite differ- 
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ent approaches. While the lSerra et al.l J2008h analys is was based on 
Monte Carlo Markov Chains (MCMC), Zahn et all us ed an approx- 
imate analytic computation of the biases (see also iHuterer et al.l 
120061) . 

In this study, we present an analytical derivation of the bi- 
ases on the cosmological parameters that goes beyond the afore- 
mentioned approximations. It is a method valid when the primary 
signal and secondary contribution (astrophysical or systematic) are 
additive and it is exact as it can be applied even when the sec- 
ondary signal is dominant over the primary. The method presented 
here is applied to the estimate of the cosmological parameters 
with the CMB power spectrum. We furthermore focus on one sin- 
gle source of additional anisotropies: those associated with unde- 
tected chisteni/nieSim (SZ) effect of galaxy clus- 
ters dSunvaev &Zel'dovicrJll972h is indeed the major source of 
secon dary temperature anisotropies ( lAghanim. Maiumdar & Silkl 
l2008t and references therein). The SZ effect is two-fold: the 
thermal SZ due to the inverse Compton scatte ring of photons 
off the hot electrons in the intra -cluster gas (e.g. iRephaelil 1 1 9951 : 
lltoh. Kohvama & Nozawa I Il998h . and the kinetic effect due to 
the Doppler shift caused by the proper motion of the clusters in 
the CMB reference frame. The upcoming large multi-frequency 
surveys will be able to detect and extract galaxy clusters using 
their specific thermal SZ spectral signature. We nevertheless ex- 
pect some level of residual SZ contribution from undetected clus- 
ters in the temperature anisotropy maps. Such a residual signal 
might be the cause of the excess of power measured by small scale 
CMB experiments like BIM A (Berkeley Illinois Maryland Associa- 
tion Array), CBI ACB AR dDawson et al.l2002l ; |Pearsori2003l ; lKuol 
l2004lReadheaj2004l) . The excess could be also due to unremoved 
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point sources dToffolatti et"ai]|2005l ; lDouspis. Aghanim & Langer 
l2006h . 

Our article is organised as follows, we present in Section 2 
the method to calculate the biases on cosmological parameters. We 
then apply our method to estimate the biases induced by the SZ 
residuals. We present our results in Section 3 and discuss them in 
the following section. Finally, we conclude in Section 5. Through- 
out the article, we assume a flat A cold dark matter (ACDM) cos- 
mological model and use for the reference model the Wilkinson 
Microwave Anis otropy Probe5 (WMAP5) cosmological parameters 
(Dunklev 2009), unless otherwise stated. 



2 DERIVATION OF THE BIASES 

In order to forecast the capability of instruments in terms of param- 
eter estimation, we commonly use the Fisher approach as it pro- 
vides us with the achievable precision on the parameters as a func- 
tion of the instrumental characteristics (noise, beam, etc.) However, 
predicting the precision is often not enough, and one should, in ad- 
dition, estimate how sensitive the measurements are to any system- 
atic or additional contribution to the signal. The aim of our study 
is to propose an analytical calculation of the bias introduced on 
the estimated parameters by the presence of an additional, non pri- 
mary signal. We consider the latter to be any signal that is added, or 
subtracted, to form the observed signal used for parameter estima- 
tion. The method presented here is thus general and can be applied 
not only to astrophysical contaminating signals (e.g. foregrounds) 
but also to instrumental systematics when they are additive. We 
develop our method in the framework of cosmological parameter 
estimation using CMB angular power spectra and we apply it to 
the case of future Planck observations in Section|3] 

Let us assume that the signal (Cf ) used to estimate cosmolog- 
ical parameters is the sum of a primary signal [Cjr MB (§)] that de- 
pends on cosmological parameters 6 and an additional signal (C'f 6 ) 
that may, or may not, contain cosmological information. If such an 
additional signal is taken into account in the parameter estimation 
analysis (Case 1) then one recovers the "true" cosmological param- 
eters, \0. If, on the contrary, one uses the total signal but assumes 
that it is made of primary only (Case 2), one obtains a set of bi- 
ased cosmological parameters, -ft. In the following, we derive an 
analytical expression of this bias, b = 2O - 16. 

We consider that the errors associated with the "data" Cf are 
distributed according to a Gaussian law. The likelihood function is 
then written as _£ = e~ x 11 . For Case 1, the;f? is : 



0) = J] cov^(CfCj) [c° x - Cf™ d '(0)] x [C° y - Cj w (8)] 

C,X,Y 

(1) 

where Cf mod ' = Cf CMB + Cf dd . For Case 2, the X \ is given by: 
X 2 2 (6) = ^ COV7 1 (C?C*) [C™ - C?™ d (0)] x [Cf - cj^m] 

CXY 

(2) 

where Cf mod = Cf CMB only, and X, Y = TT, EE, TE. We consider 
the temperature and E-mode auto-correlations, TT and EE, and the 
cross-correlation, TE, for which the coefficients of the cova riance 
matrix are cov f (C*Cj) (see, e.g. IZaldarriaga & Seliakl 199% . 

Assuming the "data" represent the sum of primordial and ad- 
ditional signals, (cf) = Cf MB ( A G) + Cf d , the set of parameters { 



minimises (x\) and the parameter set 26 minimises ix\): 
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An ensemble average of a second order approximation of equation 
l(2j gives: 
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Using equations ((3) and (0, we obtain: 
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that can be written in the compact form 
V = Gb 

wherefrom the bias vector follows simply, 
b = G 1 V. 



(4) 



(5) 



(6) 



(7) 



From equation {2), in the context of CMB angular power spectra 
measurements, 
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Equation thus allows us to calculate the bias on a parameter 
as a function of the additional signal and of the first and second 
derivatives of the primary signal. The computational advantage is 
that any additive contribution can be readily inserted without the 
need of a re-computation of the derivatives of the primary signal. 

In the case of CMB (temperature and polarisation) observa- 
tions, when an additional signal is ignored, the biases induced on 
cosmological parameters are obtained from equations (|7}{9]l. The 
derived formula (equation[5j can, however, be used in many other 
cases (matter power spectrum, weak lensing power spectrum, etc.), 
as it accounts quite generally for additive contributions or system- 
atics. Furthermore, it is the most general expression for the biases 
(for additive signals) as it applies to cases where the secondary 
signal is dominant over the primary. In co ntrast, previous studies 
(e.g. IZahn e t al. 2005? iHuterer et alj |2006) of the bias on cosmo- 
logical parameters used approximations applicable only when the 
additional signal is negligible with respect to the primary signal. 

The biases on the investigated parameters become relevant 
only if they are larger than the expected confidence intervals. The 
latter can be computed through a Fisher matrix analysis. The 68.3 
per cent confidence interval on one parameter (the others being 
known) is given by jPress et ail 1992b : 



69, = 



(10) 
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where the matrix coefficients are: 



Fa = 



dCf dCj 



Z-iZ-i c w ' dOt d6j 



(11) 



t X,Y 

The numerical values for the instrumental noise and the beam used 
to compute the covariance matrix can be found in the The Scien- 
tific Programme of Planck (also known as the Planck Blue Book, 



The Planck Collaboration 2006). 



3 COMPUTATION OF THE THERMAL SZ RESIDUAL 

The Planck satellite will measure CMB anisotropics with an un- 
precedented precision, in temperature and polarisation, over the 
full sky and from the largest scales down to five arcminutes. Fore- 
grounds and secondary anisotropics are expected to contribute to 
the signal. Taking advantage of the multi-frequency observation, 
cleaning algorithms are developed to disentangle the primary sig- 
nal from contaminants. On small scales, the SZ effect from galaxy 
clusters will be one of the major secondary contribution to the sig- 
nal. Fortunately, the characteristic spectral signature of the thermal 
SZ effect makes it easily detectable. On the one hand, one can re- 
move some of the thermal SZ signal from the CMB maps in order 
to recover the best primary CMB angular power spect rum, and on 
the o ther hand a cluster catalogue can be built (e.g. ISchafer et al.l 
2006). Nevertheless, some level of residual, unresolved SZ signal 
will remain in the temperature maps. We thus apply the method de- 
scribed in the previous section to estimate the bias induced by that 
SZ residual signal in the CMB maps on the six 'standard' cosmo- 
logical parameters. 

We now compute the additional contribution to the CMB sig- 
nal that enters equation Q. Here, it is simply the power spectrum 
of the residual SZ signal that remains after cluster extraction. We 
focus on the P oisson contribution to the SZ angular power spec- 
trum following Komatsu & Seljak ( 2002|) and assume that the con- 
tribution from correlated h alos is negligible (valid for I > 300, 
iKomatsu & Kitava ma 1999): 
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where is the comoving volume per unit redshift and solid an- 
gle and n(M,z)dM-j^ is the probability of having a galaxy clus- 
ter of mass M at a redshi ft z in the direc tion dSi. In the present 
study, we use the ISheth & Tormenl d 19991) mass function n(M,z)- 
The SZ frequency dependence is encoded in the function f(x) = 
JCpcr - 4 where x = hf\v/k^,T c is the dimensionless frequency, 
and y f = yi(M, z) is the two-dimensional Fourier transform on the 
sphere of the 3D radial profile of the Compton y-parameter of indi- 
vidual clusters, 
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with y3o(r) - erj B m °^ n e (r) and Da = Da(z) the proper 
angular-diameter distance. To model the SZ signal from individ- 
ual clusters, we assume that the electron density n c (r) within the 
cluster virial radius r v j r follows a spherical isothermal /3-profile 
dCavaliere & Fusc<> Femiand [19761) with core radius r c = 0.1 r vil . 
and B = 2/3 for all clusters for simplicity. The isothermal tempera- 
ture of the intra-cluster gas T c is set equal to the virial temperature. 

In order to evaluate equation d!2t . we choose z raax = 7 and 
A/min - 10 13 M Q . The angular power spectrum of the residual SZ is 
the contribution from undetected clusters. It is computed by setting 



A^max = Mn m (z), the cluster detection limit. The selection function 
determines the limiting mass Mn m (z) as a function of redshift for 
a cluster to be detected by a given instrument characterised by its 
beam, sensitivity and fre quency coverage. W e compute the selec- 
tion function similarly to iBartelmanrj ( l200ll) . A galaxy cluster is 
detected if at the same time its beam-convolved Compton param- 
eter y{6) emerges from the confusion noise, and if its integrated 
signal is above the instrument sensitivity. This translates into the 
following condition on the flux variation: 



AF(x) = g(x) I %G)da > A g(x) IoKii, 



(14) 



where g(x) = and Io 



2 (*b7cmb) 3 _ The integral is calculated 
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over the lines of sight for which y(8) > 3 Aybg, and the sensitivity 
limit ?i im is derived from the antenna temperature sensitivity of the 
instrument. 

The background SZ signal responsible for the confusion is 
computed assuming a Poisson distribution of clusters with masses 
between 10 13 M o and 5 x 10 l6 M o : 



Ay bK 



J <k-£jfif dMn(M,z)(^J d£ly(6,M,z)j ■ (15) 

Finally, the flux in equation <14t is integrated over a top hat func- 
tion to account for the frequency response of the instrument, and 
the integer A is the detection threshold in terms of the instrumental 
noise cry. 

In order to compute the SZ residual temperature power spec- 
trum that enters equations Q and l[8}, we derive an idealised cluster 
selection function for Planck (not accounting for the effects of fore- 
grounds and point sources) considering only the channels relevant 
for the SZ measurement, namely 100, 143 and 353 GHz (217 GHz 
is the frequency for which the thermal SZ is null). We consider 
that an SZ cluster is detected only if the selection criterion d!4t is 
satisfied in all three frequency channels simultaneously. 



4 RESULTS 

We now illustrate the application of our formalism by tak- 
ing the future Planck satellite as an example of high sensitiv- 
ity high resolution CMB experiment. The instrumental charac- 
teristics that we used can be found in the Planck Blue Book 
(The Planck Collaboration 2006). We quantify simultaneously the 
biases induced on the cosmological parameters fl A , Sl^, Hq, n s , 
cr g and r by the residual SZ contribution at 100 GHz. As a refer- 
ence model, we take the cos mological parameters obtained by the 
WMAP team tounklevl[2009T) and the corresponding C™ B + C^ es , 
where C^ cs is the residual SZ signal. 

Undetected clusters contribute as a residual signal at high mul- 
tipoles where they dominate over the primordial CMB for I higher 
than 2000 approximately (Figure Q] left panel). At I > 1000, we 
show on the right panel of figure[T]that the residual SZ contribution 
after extracting clusters above 3o> already represents more than 10 
per cent of the total signal and exceeds the instrumental noise. 

In equations l|7) and ([SJ, we do not consider polarised resid- 
uals at high £s since the polarisation induced by galaxy clusters 
is negligibl e compared to the primary CM B polarisation at those 
scales (e.g. lLiu. da Silva & Aghanim 1 12005). Therefore 



Cf e8 (0) 



Cj Tres (0) ifX = TT 
otherwise. 



(16) 



We show the induced biases in figure[2] The black solid lines 
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Figure 1. Left-hand panel: At 100 GHz, primary CMB (solid line), SZ angular power spectrum of the whole cluster population (long-dashed green line), 
residual SZ spectrum after the extraction of clusters above 3<xy simultaneously at 100, 143 and 353 GHz (dot-dashed red line), primordial CMB + residual 
SZ spectrum (dashed blue line). The black dotted lines represent the l<x error bars for the 100 GHz Planck channel. Right-hand panel: Contribution of the 
residual SZ power to the total signal after 3cry cluster detections (solid black line), Icry (green dotted line) and Scry (dashed red line). The envelop represents 
the instrumental sensitivity and the cosmic variance. 
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Figure 2. 68.3 per cent joint confidence regions for Cl\, Qj,, Hq, n s , erg and t (solid line) obtained with the expected Planck TT, TE and EE spectra computed 
for the reference cosmological model. The dotted (red) and long dashed (green) shifted ellipses represent the 68.3 per cent joint confidence regions around the 
biased parameters if respectively the icry and 5<xy residual SZ signal is not taken into account. 



represent the 68.3 per cent joint confidence level ellipses on the 
parameters Q.^, flb, Hq, n s , erg and t using the TT, TE and EE 
power spectra centered on the reference model. As shown in the 
figure, the expected constraints on cosmological parameters are at 
a few percent precision. We show, in addition, the biased values of 



the parameters due to the residual SZ signal. We consider differ- 
ent cluster detection limits (1,3, and 5cry) and thus obtain various 
levels of residuals in the signal. The results for all parameters are 
summarised in Table Q] In figure [2] the dotted red and long dashed 
green ellipses, shifted with respect to those in solid black, represent 
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Figure 3. Solid line (reference model) : Primary CMB + SZ residuals calcu- 
lated using the reference model. Dot-dashed line : Primary CMB calculated 
using the reference model. Blue dashed line : Primary CMB calculated with 
the biased parameters. The dotted lines are the envelopes representing the 
instrumental sensitivity around the reference model (CMB+residual SZ). 



Figure 4. Biases, in units of the expected lcr accuracy, induced by the SZ 
residuals as a function of the SZ cluster detection threshold for the 6 pa- 
rameters case. Aa = blue pluses, Cl b = blue stars, Hq = green x, n s = red 
diamonds, erg = red triangles and r = green squares. 



the 68.3 per cent joint confidence regions if the SZ residual signal 
that is obtained when clusters are extracted above 1 or 5a> is not 
modelled properly. 

First, and unsurprisingly, the biases induced on C1 A and H 
are negligible. At most, i.e. for a large residual contribution, they 
reach roughly 0.1 and 0.6 in units of the lcr errors on the param- 
eters. Then, as expected, cr 8 , and to a larger extent n s , fl b and r, 
are on the contrary significantly affected by the residual SZ signal. 
Those parameters are indeed, by nature, the most sensitive to an ex- 
cess of power at small scales; additionally they are degenerate. An 
excess of power at high £, due to SZ residuals, can be accounted 
for by over-estimating the value of cr 8 . Since the amplitude of the 
CMB power spectrum strongly depends on cr 8 , a relatively small 
variation of <x 8 is enough to fit the power excess. As a result, the 
bias on this parameter is rather small. The excess of power at small 
scales decreases the ratio between the amplitudes of the forth and 
fifth CMB peaks and slightly shifts them towards higher £ values 
(Fig.0 mimicking the effects induced by an increase of Q.\,. We 
find that the bias is roughly 1.4 and 2.4 in units of the error on 
Q b , for the 1 and 5cr Y cases respectively. Higher Q b and cr 8 values 
imply a higher amplitude of the primary CMB spectrum at all mul- 
tipoles. However, the SZ residual signal contributes significantly 
to the CMB spectrum only at multipoles higher than 500, leaving 
lower multipoles almost unaffected. The effects at low multipoles 
of simultaneously higher Cl b and cr 8 , namely an overall rise of the 
power, are compensated by an increase of the spectral index n s , 
which raises (lowers) the power at multipoles larger (smaller) than 
~ 1200. They are also compensated by a decrease of the optical 
depth t, that reduces the CMB power. As a consequence, fitting 
data containing primary CMB and SZ secondary residual with a 
pure primary CMB spectrum induces quite an important bias on n s 
and r. The bias on r is between 2.8 and 6 times the expected pre- 
cision for cluster detection limits 1 and 5o"y, respectively, and the 
bias on f? s is between 6.5 and 10.4 times the expected precision for 
1 and 5<Ty cluster thresholds. All these biases are summarized in 
Fig.|4]as functions of the SZ cluster detection threshold. 



5 DISCUSSION 

We show that an unremoved thermal SZ contribution in future high 
sensitivity high resolution CMB experiments, measuring both tem- 
perature and polarisation, introduces excess of power at small an- 
gular scales. This in turn induces biases on the cosmological pa- 
rameters associated with inflation, the spectral index n s and, to a 
lower extent, the normalisation erg, as well as on the density of 
baryons £1\, and the optical depth t. The other two parameters of 
the standard model remain essentially unaffected. 

The amplitude of the biases strongly depends on the ampli- 
tude and the detailed shape of the residual power spectrum. In or- 
der to estimate the residual thermal SZ power spectrum, we used 
a theoretical SZ cluster selection function built from the expected 
capabilities of Planck and not accounting for foregrounds and point 
sources. The detected clusters are expected to have masses of a few 
10 14 M G at redshifts 0.4 < z < 0.7. Although, the theoretical se- 
lection function does not coincide totally with that obtained from 
simulations (false detections, confusion with point sources, etc., J.- 
B. Melin, private communication), it gives us a reasonable estimate 
of the residual thermal SZ power spectrum used in our analysis. 

It is important to bear in mind that other astrophysical contri- 
butions need to be taken into account in future CMB experiments 
observing at small scales. In that respect, we have also estimated 
the biases induced on the cosmological parameters by the kinetic 
SZ anisotropies of the whole cluster population. Such a signal can- 
not be separated out from the primary signal as they both share 
the same spectral signature. The kinetic SZ signal, although one 
order of magnitude smaller, will dominate over the thermal SZ at 
217 GHz where the latter vanishes. We found that the bias on the 
spectral index n s is less than 2cr, whereas the other parameters are 
essentially unaltered by the additional contribution. The channel 
at 217 GHz could seem more appropriate for the determination of 
cosmological parameters. However, one needs to consider addition- 
ally the contribution from extra-galactic point s ources, radio or IR 
galaxies and their clusterin g. A recent study bvlSerra et alj (2008) 
has extended the work of iDouspis. Aghanim & Langer I ( 120061) by 
focusing on IR galaxy clustering with WMAP and ACBAR data, as 
well as with Planck. They showed (at 143 GHz) that the expected 
biases on £l\>, n s and erg are of the order of 2. 

We have also applied our method to estimate possible biases 
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Table 1. Biases induced on cosmological parameters when the SZ signal due to undetected clusters is ignored. 



detection 


detection threshold 
fiim (arcmin 2 ) 


biases expressed 
Ha fib 


in units of the expected precision 
H « s o-g r 


l<Ty 


1.6 x 1CT 4 


-0.078 


1.36 


0.3 


6.49 


0.92 -2.78 


3<Xy 


4.9 x 1CT 4 


-0.073 


2.02 


0.49 


9.12 


0.89 -4.81 


5<ty 


8.2 x 1(T 4 


-0.07 


2.36 


0.6 


10.35 


0.78 -5.96 



on the cosmological parameters with WMAP. Due to its frequency 
coverage (v < 94 GHz) we can safely consider that the thermal SZ 
residual signal in the data is constituted of the contribution from all 
the clusters (none of them being actually detected and extracted). 
The angular resolution of WMAP is limited to I < 1000. In this 
range, the contribution of the SZ signal to the CMB power does 
not exceed 10 per cent. As a result, although the contamination 
by thermal SZ is maximum in the WMAP case, we found that the 
derived cosmological parameters are unbiased. 

Given the possibility that large biases could arise due to vari- 
ous systematics, it is important to find a way to reduce their impact. 
First, one could consider only the multipoles for which the primary 
contribution to th e TT signal is dom inant, i.e. typically £ < 1000. 
This was done bv lZahn et al] (2005) who suggested to take into ac- 
count the full polarisation power spectra but to use the TT power 
spectrum only out to C ~ 1000. We have tested such an approach 
on our issue: the SZ residuals. We find that truncating the mea- 
surements above E = 1000 artificially increases the error-bars on 
Q A > f^b and Hq by approximately 10 per cent. This effect is even 
more important for fi s as it increases the error-bars by more than 30 
per cent. As expected, the cosmological parameters are then less 
biased. Nevertheless, n s remains significantly affected as the bias 
is still more than 3 times its expected accuracy. The bias on <x 8 is 
also significant. Consequently, a more appropriate analysis would 
rather consist in fitting the CMB data with a coherent model ac- 
counting for both the primary and the residual signal with full de- 
pendency on the cosmological parameters. For that purpose, one 
would need to use the cluster selection function, the understanding 
of which is rather complex (instrumental effects, foregrounds, lim- 
its due to the component separation techniques and cluster extrac- 
tion methods, etc). In that case, a better solution would be rather to 
use the total signal including the primary CMB and the secondary 
anisotropics to determine the cosmological parameters in a coher- 
ent way, that is taking into account the full cosmological depen- 
dence of the secondary signal. The la tter is crucial, as it was shown 
(Douspis, Aghanim & Langer 2006) that modelling the SZ contri- 
bution by a fixed-shape power spectrum with a varying amplitude 
Asz is not sufficient and biases the cosmological parameters. 



6 CONCLUSIONS 

In this study, we develop an analytical method to calculate the bi- 
ases on the cosmological parameters. Our method applies to any 
contamination provided the primary signal and the contamination 
are additive. Additionally, it is an exact derivation to the second 
order with the advantage of being applicable even when the con- 
taminant dominates over the primary signal. The next generation of 
CMB experiments will measure, with a high sensitivity, the signal 
at small angular scales where secondary contributions intervene. 
We apply our method to the case of a contribution from undetected 



thermal SZ clusters to the primary CMB signal (assuming no con- 
tribution from foregrounds or point sources). 

For illustration, we take the characteristics (noise, beam size) 
of Planck and compute the residual SZ signal from the undetected 
clusters assuming all clusters above 3 or 5cr Y are detected simul- 
taneously in the channels 100, 143 and 353 GHz. The residual SZ 
signal contributes more than 10 per cent of the total signal at multi- 
poles higher than 1000. The higher the SZ cluster detection thresh- 
old, the higher the contamination. 

We perform a bias estimation simultaneously on six cosmo- 
logical parameters (Q.a, Clb,Ho, n s , o~& and r) using temperature and 
polarisation anisotropy TT, TE and EE power spectra. This quanti- 
fies the effect of fitting the data, that include a residual contribution, 
with a model that ignores it. We then compare the biases to the ex- 
pected lcr errors on each parameter. We find that the biases induced 
by the thermal SZ residual signal on Q.a and Hq are negligible. At 
most, they are of the order of 0.08 and 0.6 in units of the lcr error 
on the parameters. On the contrary, the determination of Q. b , n s and 
t is significantly altered by the residual SZ signal. The biases are 
2.4, 6. and 10.4cr respectively. This is easily understood as they are 
the most sensitive parameters to an excess of power at small scales 
and, moreover, they are degenerate. 

We point out the importance of taking into account the SZ 
residuals in the analysis of the small scale high sensitivity CMB 
data. The SZ residual power spectrum depends on the cosmological 
parameters and on the cluster selection function. A joint analysis 
of primary and secondary CMB signal will provide additional con- 
straints on the cosmological parameters and thus reduce the biases 
arising from the SZ residual power excess at high multipoles. A 
coherent analysis, including full cosmological parameters depen- 
dency of the primary and secondary signal, allow one to use the 
whole range of multipoles, including the highest ones. 
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